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We study the relation between center vortices and instantons in lattice QCD. 



1. Motivation 

One approach to understanding the physical 
origin of confinement and chiral symmetry break- 
ing is to attempt to capture the essential physics 
in a much smaller number of degrees of freedom 
by projecting from a maximal abelian or maximal 
center gauge. For simplicity, we consider SU{2) 
and denote a link variable transformed by a gauge 
transformation g by 



U^x)=g\x)U^x)g{x + l j) 



= („9 



k 



(1) 



with — 1 ■ Maximal abelian gauge corre- 

sponds to calculating 

Max {g} ^tr(l^)a 3 [/fOr)a 3 ) 

= Max {s} ^{2[(ag) 2 + (af) 2 ] S)M -l} 

and the abelian projection is performed by trun- 
cating ai and ai and normalizing to obtain a U(l) 
matrixEJ. Although results are ambiguous be- 
cause of Gribov ambiguities associated with local 
minima, the projected results generally reproduce 
the proper string tension. 

Maximal center gauge corresponds to calculat- 
ing 

Max {9} ^(TrC/9( x )) 2 = Max {g} ]T 4(a 2 )^ 
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and center projection is performed by truncat- 
ing ai, «2i an( i a 3 an d normalizing the resulting 
Z(2) variable to ±1. The properties of center pro- 
jection have been discussed widely by Greensite 
and collaborators [|j and Reinhardt and collabo- 
rators^]. The set of negative links defines a set 
of center vortices, corresponding to closed lines 
in three dimensions and closed surfaces in four 
dimensions. Vortices may also be thought of as 
Dirac strings connecting monopoles. 

For subsequent reference, it is useful to visu- 
alize a vortex sheet generated by a single neg- 
ative link in the z direction from (0, 0, 0, 0) to 
(0,0, 1,0). The vortex sheet is the surface of the 
cube dual to this link, that is, the unit cube in 
(x,y,t) centered at the origin and having z=\. 
If one cuts this four-dimensional surface at time 
t = and observes it in the three-dimensional 
x, y, z space, it corresponds to a square loop pierc- 
ing the center of each of the plaquettes with Wil- 
son loop —1 that share the negative- z link. If two 
adjacent links in the z direction are —1, the cube 
grows to a rectangular parallelepiped, and in gen- 
eral surfaces can grow by stacking the dual cubes 
together. Returning to our single cube dual to 
a z link, we note for example that the x-t sur- 
face corresponds to a y-z plaquette and thus cor- 
responds to B x ~ F yz . In general, the vortex 
sheet in the ij plane corresponds to SijkeFke- By 
suitably choosing the time direction in euclidean 
space to either be in or perpendicular to the plane 
of the sheet, we can think of the vortex sheets as 
being sheets of magnetic or electric fields. Hence, 
a single smoothly bending vortex sheet corre- 
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sponds locally to either E or B, but can never 
generate E ■ B. The simplest way to generate 
lumps of E ■ B, corresponding to the topological 
charge density one expects from instantons, is to 
have an intersection of two mutually perpendicu- 
lar sheets. For example, a sheet in the x-y plane 
at fixed (to, zq) with F z t ~ E z intersecting a sheet 
in the z-t plane at fixed (xo,yo) with F xy ~ B z 
yields a point of E ■ B at (xq, yo, z , t ). If the 
vortices have finite thickness, the resulting lump 
of E ■ B will have finite size. 

Lattice phenomenology provides several in- 
dications that the reduced degrees of freedom 
subsequent to center projection still embody 
the essence of confinement and chiral symme- 
try breaking. The contribution of each vor- 
tex piercing a Wilson loop is (—1), so that a 
random ensemble of vortices produces confine- 
ment and lattice measurements yield the proper 
string tension and scaling ^|J3). There are sev- 
eral indications that there are underlying physical 
structures ( "thick vortices" ) corresponding to the 
"thin vortices" obtained by projection. Calcula- 
tions with twisted boundary conditions to pre- 
vent decay to the vacuum yield stationary vortex 
solutions that are exponentially localized in two 
dimensions and periodic in the other two dimen- 
sions Q. When the contribution of a lattice gluon 
configuration to a Wilson loop W n is labeled by 
the number n of projected vortices that would 
pierce the loop, one observes that [W2)/(Wq) — » 1 
and (Wi)/(W ) -> -1 for large loops [g. This 
suggests that there is a finite size object in the ac- 
tual configuration that, when enclosed by a large 
enough loop, behaves just like an idealized "thin 
vortex" . In addition, at finite temperature, the 
vortices appear to become aligned in the time di- 
rection, so that they pierce spatial loops yielding 
an area law. Finally, de Forcrand and D'Elia|| 
have studied the effect of removing center vor- 
tices. For each link U^x) one determines Z^[x) 
by center projection and defines the modified dis- 
tribution U'(x) — Z fJ _(x)U fJ _(x), which projects to 
the trivial configuration. Not only does one lose 
confinement, but in addition (ipip) — > and the 
topological charge (Q) — > 0. That (Q) — > is un- 
derstandable, since U must be far from 1 some- 



where for nontrivial topology and projection leads 
to the trivial sector. That — * implies via 
the Banks Cacher relation that the zero modes 
and thus all the instantons and antiinstantons 
must also be removed. 

2. Center Projection of a Single Instanton 

To understand how the modified configuration 
removes instantons, we have studied the center 
projection of a single instanton on a lattice. Pro- 
jection yields a cluster of parallel —1 links in the 
center of the instanton in one common direction 
that we will call time. This breaks the euclidean 
invariance of the instanton, and yields a vortex 
sheet that is pure E. The most relevant config- 
urations are a single dual cube, a 2 3 cube com- 
posed of 8 elementary cubes, and a 3 3 cube com- 
posed of 27 elementary cubes. To see how the 
modified distribution cancels an instanton, note 
that there must be a shift in group space of ir so 
U ~ e iE -> e i( - E+n l Since 
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if we assume that the full shift of it occurs in 
a single plaquette for an instanton of size p and 
assume a spherical sheet of radius i?, 
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so that 5Q — 1 when R = p. For a cube, SQ = 
1.203 for R = 0.813 p. Hence, if the size of the 
vortex is ~ p, it should cancel Q. 

The primary technical problem is local minima 
corresponding to Gribov ambiguities. We start 
with a cooled instanton on lattices from 16 4 to 
32 4 , minimize 

# {9} = ^[1- (±tTg\x)U^x)g{x + ii)f] 

successively on each site by matrix diagonaliza- 
tion, and overrelax by applying g a with a *~ 1.7 
to 1.8. The local minimum depends upon a. For 
example, a regular gauge instanton with p = 4 on 
a 16 4 lattice yields no vortex for a = 1.7 and 
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Figure 1. Spherical average of topological charge 
density times r 3 versus radius r. The solid curve 
denotes a cooled instanton, the crosses indicate 
the density after removal of a 4 3 vortex, and the 
squares and stars denote the density after one and 
two cooling steps, respectively. 



cancel the original toplogical charge, the vortex 
sheet by itself has no topological charge and thus 
fails to represent the essential phyiscs of the orig- 
inal configuration. 

As argued above we expect lumps of E ■ B to 
arise from intersections of two vortices. Whereas 
an isolated instanton will not have vortex sheets 
extending to infinity, it would be natural for vor- 
tices to connect neighboring instantons and anti- 
instantons. Thus, the relevant case to examine is 
an ensemble of instantons and antiinstantons, and 
to see whether projection produces intersections 
with localized instanton- like lumps of E ■ B. Were 
this to occur, center vortices would succinctly rep- 
resent the essential physics of both confinement 
and chiral symmetry breaking. 
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a 2 3 vortex for a — 1.8. Similarly the initial 
gauge matters. A random gauge transformation 
at a = 1.7 yields a 2 3 vortex whereas the regular 
gauge does not. Also, going first to Lorentz gauge 
can prevent finding a vortex solution. By explor- 
ing a variety of gauges and initial conditions, the 
2 3 vortex emerged as the minimum solution for 
this case. 

As expected from the analytic estimate of 
(5Q Lat , for a vortex size R comparable to p, cen- 
ter projection removes the instantons. Figure [l] 
shows r 3 times the topological charge density as 
a function of distance for a cooled instanton of 
size p — 4.7, the instanton with a vortex of size 
R = 4 removed, and the density after one and 
two cooling sweeps. Since 6Q begins localized on 
a shell, it takes several cooling sweeps to see that 
the integral averages to zero. After more cooling 
sweeps, the entire density relaxes to zero. A sim- 
ilar calculation with size R — 2 cools back to the 
original instanton with Q = 1. 

With hindsight, we believe this single- instanton 
example is an unrepresentative special case. The 
problem is that the vortex sheet is a pure E field. 
Although when it has the proper area, it can com- 
bine with the existing B field of the instanton to 
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